Preserving qubit coherence by dynamical decoupling 
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In quantum information processing, it is vital to protect the coherence of qubits in noisy environments. 
Dynamical decoupling (DD), which applies a sequence of flips on qubits and averages the qubit-environment 
coupling to zero, is a promising strategy compatible with other desired functionalities such as quantum gates. 
Here we review the recent progresses in theories of dynamical decoupling and experimental demonstrations. 
We give both semiclassical and quantum descriptions of the qubit decoherence due to coupling to noisy 
environments. Based on the quantum picture, a geometrical interpretation of DD is presented. The periodic 
Carr-Purcell-Meiboom-Gill DD and the concatenated DD are reviewed, followed by a detailed exploration of 
the recently developed Uhrig DD, which employs the least number of pulses in an unequally spaced sequence 
to suppress the qubit-environment coupling to a given order of the evolution time. Some new developments and 
perspectives are also discussed. 
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I. COHERENCE AND DECOHERENCE 

The power of quantum information processing y]], the 
quantum parallelism, comes from the superposition principle 
of quantum mechanics. The building block of quantum tech- 
nology, a quantum bit (qubit), is a two-level system that can 
be identified as a spin-1/2 with states |f) and |X). The ability 
of the qubit to be in a coherent superposition of |f) and |X), 



l^f) = cos —e 



^ 2 |T> + sin- e '^ 2 |i), 



(1) 



enables the parallel processing of many pieces of classical in- 
formation. In order for this idea to work, the qubit has to 
faithfully maintain its quantum state. Not only the populations 
cos 2 (6>/2) and sin 2 (6»/2) in the two states |f> and ||>, but also 
the relative phase e~ ,lf between |T) and |X) should be kept at 
certain values. Unavoidable couplings between the qubit and 
the environment (hereafter referred to as bath) spoil the quan- 
tum state by introducing uncontrolled evolution of the qubit. 
The populations and phases are randomized and the qubit co- 
herence is lost. This decoherence problem is one of the most 
serious obstacles in the roads towards scalable quantum infor- 
mation processing J2]. 

The population randomization (i.e., relaxation) process in- 
volves energy dissipation and therefore is subjected to the en- 
ergy conservation condition. Thus it can be suppressed by in- 
creasing the spin splitting of the qubit. In contrast, the phase 
randomization (i.e., pure dephasing) is a more serious issue, 
since this process does not involve energy dissipation. 



A. Semiclassical picture of decoherence 

In the semiclassical picture, the pure dephasing of a qubit 
or a spin-1/2 is caused by the fluctuation of a local classical 
field fixed at a given direction OlH]. The Hamiltonian of the 
qubit in the external field including the random component is 



H, 



qubit 



[wo/2 + Z(0] & t , 



(2) 



where & z is the Pauli matrix for the qubit, «o is the Zeeman 
splitting under the external field, and 2Z(t) is the random field 
resulting from the interaction with the bath. Let us consider a 
qubit initially in a coherent superposition state 



|^(0)> = C + |T> + C_IX>, 

corresponding to a pure state density matrix 



P(0) 



\c + 



c+c*_ 



C* + C- \C- 



(3) 



(4) 



in the basis |f ) , |X). At the end of the evolution, a random 
relative phase <p(r) — 2 jC Z(t)dt between |f) and |X) is accu- 
mulated in the qubit wave function 



[<A(r)> = C + e 



-iip(T)l2 



IT) + C_ e 



,y.(T)/2 



II): 



(5) 



and the off-diagonal coherence of the resulting density matrix 



P(t) 



IC+I 



C+C . e -,V(T) 

\C-\ 2 



(6) 
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becomes random. The ensemble average over all possible 
realizations of the random noise Z(t) gives the decay of the 
off-diagonal density matrix elements, i.e., the decoherence of 
the qubit (or the depolarization of the spin-1/2 in the plane 
perpendicular to the external field). The resulting qubit state 
is a mixed state with vanishing off-diagonal coherence, since 
the noise-averaged quantity (e~" p{j A vanishes in the long time 
limit. 



B. Quantum theory of decoherence 

In the quantum picture 15[], the decoherence of a qubit re- 
sults from the qubit-bath entanglement, which is established 
during the evolution of the interacting qubit-bath system. The 
general pure dephasing Hamiltonian has the form 



H dp -C + & z 



(7) 



where C is the interaction within the bath and Z is the bath 
operator representing the quantum field on the qubit resulting 
from the qubit-bath interaction. Suppose the initial state of the 
qubit-bath system has the form |*F(0)) = [(HO)} ® \J) , i-e., a 
direct product of the qubit state \ifr(0)) = C + |f} + C ||> and 
the bath state |/). At the end of the evolution, an entangled 
state is established as 



|Y(T)> = C + |+>®e 



-i(C+Z)T 



\J) + C.\-)< 



-KC-Z) 



T \J) 



= C + |+) ® |/ (+, (r)) + C |-> ® \J ( -\t)) ., 



(8) 



and the off-diagonal coherence of the reduced density matrix 
of the qubit becomes bath-state-dependent 



P(t) = 



\C + 



C\C (7«(t)|JW(t)) 



C + C1 {j ( ->(t)\J^(t)) 



\c. 



(9) 



The off-diagonal qubit coherence is reduced when the bath 
state overlap decreases 



Lj(r) = (j ( -\t)\J (+ \t)) = (j\e h 



(c-z)T-i(c+zy 



V). (10) 



A transparent physical meaning of this formula is that the co- 
herence of the qubit decreases when the distinguishability of 
the bath states increases, or the quantumness of the qubit de- 
cays when it is gradually "measured" by the environment. 

The decoherence in Eq. ( TTUb is caused by the quantum fluc- 
tuation of the local field for a single bath state \J). At fi- 
nite temperature, the bath itself is in a thermal ensemble as 
IZj Pj\J)(J\- Ensemble average over the distribution of the 
initial bath states \J) causes additional dephasing due to the 
thermal fluctuation, referred to as inhomogeneous broadening 
in literature Ja] . 

As an example, in a confined solid-state environment such 
as a quantum dot, the most relevant source of decoherence at 
low temperature (a few Kelvin) for an electron spin is the hy- 
perfine interaction with the lattice nuclear spins (which serve 
as the bath) |g-ca]- In a moderate (> 0.1 Tesla in GaAs quan- 
tum dots) external magnetic field, the electron spin relaxation 
is strongly suppressed IU0l4l2ll and the coherence decay is 
dominated by pure dephasing. Recently, a variety of quantum 
many-body theories have been developed to evaluate the bath 
state evolution Lj{t) or its ensemble average, including the 
density matrix cluster expansion 11131, 1 1411 . the pair-correlation 
approximation J5D, the linked-cluster expansion 1131 . and the 
cluster correlation expansion Illa.[l7ll . In the pair-correlation 
approximation J5[], each pair- wise flip-flop of the nuclear spins 
is identified as an elementary excitation mode and is taken as 



independent of each other. To study the higher order corre- 
lations, the Feynman diagram linked-cluster expansion is de- 
veloped 1131 . The evaluation of higher-order linked-cluster 
expansion, however, is tedious due to the increasing number 
and complexity of diagrams with increasin g the interaction or- 
der. The density matrix cluster expansion IIO, Il4fl provides a 
simple solution to include the higher-order spin interaction ef- 
fects beyond the pair-correlation approximation (without the 
need to count or evaluate Feynman diagrams). However, the 
accuracy problem (even when the expansion converges) limits 
the cluster expansion to applications in large spin baths. The 
cluster-correlation expansion Im \Y% bears the accuracy of 
the linked-cluster expansion (the results are accurate when- 
ever converge) and the simplicity of the cluster expansion 
(without the need to count or evaluate Feynman diagrams), 
while free from the large-bath restriction of the cluster expan- 



II. SUPPRESSING DECOHERENCE BY DYNAMICAL 
DECOUPLING 

Since qubit decoherence results from uncontrolled evolu- 
tion due to the coupling between the qubit and the bath, a natu- 
ral idea to combat decoherence is to encode the qubit in a sub- 
space immune to noises from the environment (decoherence- 
free subspace IU8l 11911 ). which is made possible by symme- 
tries of the interactions in certain physical systems. Or alter- 
natively, the coherence can be protected by dynamically elim- 
inating the qubit-bath coupling during the evolution (dynami- 
cal decoupling, referred to as DD for short). The DD schemes 
were originated from the Hahn echo 112011 and were developed 
for high-precision magnetic resonance spectroscopy J2ll - l23ll . 
When the field of quantum computing was opened up, the 
idea of DD was introduced to protect qubit coherence 11241 — 
12711 . which stimulated numerous theoretical studies on exten- 
sion an d ap plications of the DD approach to quantum com- 
puting |28h38I1 . The recent experimental advances are also 
remarkable 09H41I1 . 

In the DD scheme, a sequence of pulses is applied to flip the 
qubit and average the qubit-bath coupling to zero during the 
evolution. It is a promising strategy due to its compatibility 
with other desired functionalities such as quantum gates J42l - 
l44fl . The most general Hamiltonian describing the coupling 
between a qubit and a bath reads 



H - C + & x ® X + & y ® Y + & z ® Z, 



(11) 



where d" A / v /z are the Pauli matrices for the qubit, and 
C,X, Y, and Z are bath operators. The off-diagonal coupling 
(&x ® X + & y ® Yj induces population relaxation. The diago- 
nal coupling & z ® Z induces pure dephasing. 



A. Carr-Purcell-Meiboom-Gill DD 

For the sake of simplicity, we first consider the pure de- 
phasing case (X - Y - 0). In the absence of controlling 



pulses, the evolution of the quantum state |\P(t)) = Uq ^(O)) 
of the coupled qubit-bath system is driven by the free propa- 
gator Uq = e 



iHr 



- e 



-i(C+fr z & 



The Hahn echo 12011 is realized by a single instantaneous n 
pulse applied at the middle of the evolution to switch the qubit 
states between |f ) and |X), 



PF(2t)> = U & X U mO)} , 



(12) 



so that the propagator for the whole evolution from to 2t is 

Uq& x Uq = 6~ X U\ with 



G l s & x Uq& x Uq = e -KC-&jmr e -KC*&ji 

_ -i2r[C+& z <SZ-0(Cr)] 



(13) 



In the propagator, the qubit-bath coupling is eliminated in 
the first order of the pulse interval r. By repeating the 
Hahn echo propagator U\, the Carr-Purcell-Meiboom-Gill 
DD (CPMG) B45l |46y can be constructed so as to preserve 
the coherence of the qubit for a long time. 

The building block of CPMG consists of two instantaneous 
n pulses applied at t and 3t, respectively. At the end of 
the evolution t - At, the state of the qubit-bath system is 
| x P(4t)) = U 2 |*P(0)), where the propagator 



U 2 - Uo& x UqUq& x Uq = d- x U\& x U\ 

- „-i4T[C+a- : 0ZO(C 2 T 2 )] 



(14) 



is obtained by embedding U\ into the basic structure 
d~x (' ' ' ) o".t (' ' ' )■ The CPMG sequence of 2N pulses is ob- 
tained by repeating the building block U 2 for N times. The 
propagator for the whole evolution from to T — 4Nt is 



t^CPMG - U 2 - e 



-iT[C+fr : ®Z-0(C 2 T 2 )] 



(15) 



The qubit-bath coupling is eliminated up to the second order 
of the minimum pulse interval r. 



B. Concatenated DD 

Note that in Eq. dT4j, the building unit of CPMG can be 
viewed as a nested application of the Hahn echo, which elim- 
inates the qubit-bath coupling to one order higher than the 
simple Hahn echo does. It was noticed in Ref. 12a, 12711 that 
a mirror-symmetric arrangement of two DD sequences can 
decouple a quantum object to a higher order. And further- 
more, Ref. 12711 mentioned the possibility of realizing DD 
to an arbitrary order by iterative construction. Khodjasteh 
and Lid ar propo sed the first explicit concatenated DD (CDD) 
scheme OCX 13111 to eliminate arbitrary qubit-bath coupling (in- 
cluding both diagonal and off-diagonal couplings) with an 
intuitive geometrical understanding 14711 . The idea of CDD 
was further developed by incorporation of randomness into 
the sequence for improvement of performance 11321 13611 . CDD 
schemes against pure dephasing were investigated for electron 
spin qubits in realistic solid-state systems with nuclear spins 
as baths 133113511 . The advantage of CDD over the periodic 



DD sequences has been observed in experiments for nuclear 
spin qubits in solid state environments [43|. 

The propagator for CDD is obtained by recursion 



U„ - & x U n .\& x U n . 



,-/2"T[C+o- ; i8ZO(C"r")] 



(16) 



in which the qubit-bath coupling has been eliminated up to 
the nth order of the minimum pulse interval r. By increas- 
ing the concatenation level n, the qubit-bath coupling can be 
eliminated up to an arbitrary order of r. 

For the most general qubit-bath Hamiltonian in Eq. ( fTTT ). 
the idea of concatenation can still be applied to eliminate both 
the pure dephasing term o% ® Z and the relaxation term & x ® 
X + & y ® Y . In the absence of controlling pulses, the evolution 
of the qubit-bath system is driven by the free propagator Uq = 
e -iHr -pjj e q U bit-bath coupling can be eliminated up to the 
first order of t by the controlled evolution 113011 



U\ =U [& X U & X ] \fryU &y\ [o-jf/oo-z] 



-iT(C+cr x ®X+fry®Y+fr z i 



-iT(C-& x ®X+o- y ®Y-& z < 



-i4T[C+0(a/3r)] 



-it(C+o- a -®X-o- v < 



-iT(C-d-,®X-d- y ®Y+fr z e 



(17) 



where a denotes the norm of C and /3 denotes the norm of 
X, Y, Z. Thus all the qubit-bath coupling terms are eliminated 
in the first order of the minimum pulse interval t. By concate- 
nation, the propagator for the «th order CDD is 

U„ = U n -i [o- x U„-i& x j \fr y U n -\&y\ \& z U n -i& z ] 

_ -i2"T[C+0(aP'i")+O{a 2 l3"- ] T")+--+O(a"l3T")] /ig\ 

in which the qubit-bath coupling has been eliminated up to 
the nth order of r. By increasing the concatenation level «, 
the qubit-bath coupling can be eliminated up to an arbitrary 
order. 

To eliminate the qubit-bath coupling to a given order JV 
of the evolution time, the number of instantaneous n pulses 
scales exponentially with the order of CDD, namely, A^ pu i se = 
0(2 N ) for eliminating the pure dephasing term and A^ pu i se = 
0(4") for eliminating all qubit-bath couplings. Since errors 
are inevitably introduced in each n pulse, it is desirable to 
minimize the number of controlling pulses used to achieve a 
given order of decoupling. 



C. Uhrig DD 



Uhrig DD (UDD) 137J I48U50H is a remarkable advance in 
the DD theory. UDD can eliminate the qubit-bath pure de- 
phasing up to the Mh order of the evolution time using Af 
instantaneous n pulses applied at 



T ; — Tsin" 



W 



2(N + iy 



(;=1,2,...,A0, 



(19) 



during the evolution of the qubit-bath system from to T. 
UDD is optimal in the sense that it uses the minimum num- 
ber of control pulses for a given order of decoupling. Such 



pulse sequences for N < 5 were first noticed by Dhar et al 
in designing control of the quantum Zeno effect 15111 . The 
application of such sequences to DD was first proposed by 
Uhrig for a pure dephasing spin-boson model 113711 . Then Lee, 
Witzel and Das Sarma conjectured that UDD may work for a 
general pure dephasing model with an analytical verification 
up to N — 9 14811 . Later, computer-assisted algebra was used 
to verify the conjecture up to N = 14 14911 . Finally, UDD was 
rigorously proved to be universal for any order N J5C ]and was 
also extended to the case of population relaxation 15011 . The 
performance bounds for UDD against pure dephasing were 
also established 15211 . 

By construction, UDD is optimal for a finite system (or a 
system with a hard cut-off in the spectrum) in the "high fi- 
delity" regime where a short-time expansion of the qubit-bath 
propagator converges. For the "low fidelity" regime, further 
theoretical work 08115311 shows that UDD is optimal when the 
noise spectrum has a hard cutoff, while CPMG performs better 
than CDD and UDD when the noise has a soft cutoff or when 
the hard cutoff is not reached by the spectrum filtering func- 
tions corresponding to the DD sequences. The experimental 
investigations of UDD were carried out in an array of ~ 1000 
Be + ions [40] and in irradiated malonic acid crystals 14111 . 



By exact diagonalization of the spin-boson Hamiltonian, 
6U(T) has been evaluated as 6U(T) = e lK(T) with 



ACT) = £(-iy [&i(Tj) - Kj{T j+l )] . 

7=0 

where we have defined 7q = 0, 7V+i = T, and 



Ki{t) = e iCl Ke~ iCt , 






The Taylor expansion 



dty 



/7-fold commutator 



(24) 



(25a) 
(25b) 



m) = t + J) ^-f- [c, ■ ■ ■ [c, [c, k]]] = k + Y, K P f, 
p=l I ■ p=\ 



(26) 



yields A(T) = - £ K P T''A P , where 
P =\ 



1. Spin-boson model: discovery of UDD 

The qubit-bath Hamiltonian H& of the spin-boson pure 
dephasing model 113711 corresponds to C — J^jCOthhi and 

Z = 2i( A: '/2)(^' + bj) in Eq. (0, where hi is the bosonic an- 
nihilation operator. For N instantaneous n pulses applied at 
T\, T%, ■ • ■ ,Tpf € [0, T], the propagator for the evolution from 
to T is 

U(T,0) =O (T - T N )fr x Uo(T N - T N .{) 

■■■& x U Q {T2-T l )& x lJ (T x ), (20) 

where Uo(t) = e *' is the free propagator. U{T,Q) can be 
written as Un (for N being even) or & x Un (for N being odd) 
with 

ft =e -KdH- 1 fd- : ®Z)(T-T N ) e -KC+(- 1 ) N - X & Z <SZ){T N -T N ^) 
„-i(C-fr : <8Z)(T 2 -T x ) -i(C+&-%Z)T x 



U (+) + U^ U (+) - C7 (_) 

■ h CT, ® , 



(21) 



where U^ is obtained from Un by replacing fr z by ±1. In the 
Mh order UDD, the positions T\, T%, ■ ■ ■ , Tn of the N pulses 
are fixed by requiring that in the propagator Un, the qubit-bath 
coupling should be eliminated up to the Mh order, i.e. 

6U(T) = U { + } - U^ = U^KU^'U^ - 1] ~ 0(T N+1 ), 

(22) 

or equivalently, 

6U(T) = (f/^ , ) t f> (+) - 1 ~ 0(T N+1 ). (23) 



-§-im-(?) 



(27) 



Thus the condition §U(T) - 0(T N+l ) is equivalent to N cou- 
pled algebra equations 



A^ = 0, (p=l,2, 



,N). 



(28) 



whose unique physical solution is the UDD sequence in 
Eq. ( fT9l . The UDD sequence is optimal in that it uses the 
minimum number of pulses to make the first N terms of A p 's 
vanish and eliminate the qubit-bath coupling up to the Mh 
order. 



2. Geometrical interpretation of decoherence and DD 

Here we give a geometrical interpretation of decoherence 
and DD by considering the spin-boson pure dephasing model, 
based on trajectories of bath quantum states in the Hilbert 
space conditioned on the qubit states and DD control. The 
pure dephasing qubit-bath Hamiltonian can be reformulated 
as 



// = ^|±X±|®£ ± , 



(29) 



where ||±)} denote the two eigenstates of the qubit, and the 
bath operators H ± — C + Z. The qubit coherence is given by 
the overlap of bath states, as shown in Eq. ( fTOb - 

The state of the bosonic bath can be described in the basis 



of coherent states 15411 . The coherent state of the fth boson 
mode is \Pi) = e Plb i~ p i h '\0) with P/ being a complex number. 
A coherent state |P/(?o)) after a time of evolution under the 



---.A-CTj) 



Pi,+<J 3 ) 




FIG. 1: The bifurcated trajectories (solid curves) of Pj,±(t) in the 
complex plane under qubit flip control applied at T\ = % and T 2 = 
2£, with the initial time r = and final time T } = T. The arrows on 
the solid curves indicate the evolution directions. 



Hamiltonians H/ ± = co/b^bi + jKi(b, + bf) is still a coherent 
state 



-iH L± (t-T ) 



\Pi(T )) = \Pi, ± {t))e 



'id: 



(30) 



where 



PiAt) = [Pi(T ) - (+^L)]e-^'- T ^ + £-, (31) 

and the phase factor ± (t) - ± J' %[\KiP*,{t)\dt. 

As illustrated in Fig.Q] the complex numbers Pi t ±(f), which 
represent the coherent states, are rotating clockwise about the 
points +^ in the complex plane with an amplitude of |P/(7o)| 
and an angular frequency a)/. The overlap of the bifurcated 
states 

\(Pl t+ (t)\PUt))\ = expHP/.+W - P,,_(0l 2 ), (32) 

is determined by their distance in the complex plane. We 
consider the case that the initial bath state is a coherent state 
|/) = (^). \Pi(Tq)). Thus at time t the qubit coherence 



Lj(t) = (g) \(P,At)\PiAt))l 



(33) 



decreasing when the distance between P/,+(f) and P/ _ in the 
complex plane is increased. Since the bifurcated evolution 
of the bath is determined by the qubit states |±), during the 
qubit-bath evolution, instantaneous flips of the qubit states 
will cause the bath evolution pathways to exchange their rota- 
tion centers. At some later time, the two bifurcated pathways 
could cross into each other, upon which the qubit and the 
boson mode become disentangled. At this disentanglement 
point, the which-way information is erased, and therefore the 
qubit coherence is recovered. 

Let the bifurcated bath states at time T m -\ be denoted by 
the complex numbers {PiAT m -i)}- Suppose there is a qubit 



flip applied at t — T m -\. After an interval of evolution, the 
bath states will become 



P,,±(T m ) = [P,AT m -i) ± -^k<^»-^-.> + ^L 

loj[ 2u>i 



(34) 



We define the difference A l m = Pi+iJm) — Pu-(Tm)- By recur- 
sively using the initial condition 

P,, + (T ) = PiATq) = Pt(To), (35) 

and Eq. d34b , we have that after Af flips at times T\, T2, • • • , 
TV, the difference 



A'. 



N+l 



K-D 



N+l-ica,T Ntl 



Klf{Ul), 



with 



1 N 
f(u,) = — V(-iy(e !W;r ^ -e iulT ' 

l0Ji U 



(36) 



(37) 



Eqs. (l34l and (l36l l give us a geometrical interpretation of con- 
trol of decoherence by qubit flips. In Fig. Q] we show the 
evolution of PjATm) for qubit flips occurring at T\ — j and 
T 2 = Ij-, with the total evolution time T3 -T . 

Note that the initial Pi(To) is canceled in the expression of 
A^ +1 in Eq. (l36"l l. Thus from Eqs. d32l and d33l , the coherence 
Lj is independent of the initial bath state \J) = (££), \Pi(Tq)). 
By the expansion of f(a>i), we obtain 



AL, = (-D 



N+\„-UO[T K l 
0>i 



E 



(mT) n 



A„ 



(38) 



' ! between 



where A„ is given by Eq. d27| i. The distance A 
PiATn+i) is a small quantity ~ 0(t n+1 ) if {A„ = 0} for n 
N. Thus the conditions for UDD are reproduced. 

3. Proof of universality of UDD against pure dephasing 

The proof of the universality (i.e., model indepen- 
dence) IdOII of UDD is facilitated by the observation that to 
eliminate the qubit-bath coupling to a given order, one needs 
only to eliminate the odd-power terms of the coupling <x-®Z in 
the perturbative expansion of the propagator, since the even- 
power terms of & Z ®Z is a pure bath operator, (& z ®Z) 2m - Z 2 '*, 
which does not cause qubit decoherence. We will present the 
proof in the interaction picture following Ref. 115011 . which can 
be easily reformulated in other pictures i52ll . 

As discussed in the previous subsection, for the pure de- 
phasing Hamiltonian in Eq. (0 under the control of the Mh 
order UDD sequence, the propagator from to T is given by 



(J ^ e -i(CH-r) f 'a- : 0Z)(T-T N ) e -i(CH-l) N -^ z 'SZ)(T N -T N _ l ) 
-i(C-fr : ®Z)(T 2 -T,)„-HC+& z ®Z)T 1 



(39) 



Proof of the universality of UDD is equivalent to proving 

N = ( ^ th) + 0(T N+ \ (40) 



where 0^ is a bath operator containing no qubit operators. 
With the standard perturbation theory in the interaction pic- 
ture, Eq. ( 139b can be put in the time-ordered formal expression 



N = e -£T g e -i£ FN<t)&&ZiMt^ 



(41) 



where 3f is the time-ordering operator, the modulation func- 
tion F N (t) = (-iy for t € [r ; -,r, + i] with To = and 
TV+i = T, and 



p-folds commutator 



(ity\ 



z /Ws ^z^^^[c,[c,..[c,z]...]] 



/>=<> 



p=0 



(42) 



The propagator can be expanded into Taylor series 

CO 

N = e- iCT Yji-id-zT ® A„ = ^ even) + f/J; dd) , (43) 

n=0 

where 

A„ s I F N {t„)dt„ I F N (t„-i)dt„-i 
Jo Jo 

••■ F N (ti)dt 1 2 I (t n )2 I (t„-i)---2 I (ti), (44) 
Jo 

is a pure bath operator. Here 

oo 
£(eve„) =e -,cr^ ( _ 2^ (45) 



i=0 



consists of even powers of the qubit-bath coupling & z ® Z and 
therefore is a pure bath operator, which does not induce qubit 
dephasing. The term consisting of the odd powers of the qubit 
bath coupling 

oo 

U<-° dd) = cr z ® ?** 2(-0 a+1 A 2 , +1) (46) 

fc=0 

induces the qubit dephasing. We just need to show A2t + i = 

o(t n+1 ). 

Using the expansion in Eq. ( l42b . we have 

A„ = J] %< ■ ■■Zp^P> F P<^--,pJ" +p ' +P2 --- +p "] , (47) 
[pj) 



where 



'— r*-rfrTn"W(r 



(48) 



is a dimensionless constant independent of T. Now the prob- 
lem is reduced to proving 



Fpi,P%,-,Pn — ^' 



(49) 



for n being odd and n + Yj%\ Pj ^ N- F° r this purpose, we 
make the variable substitution tj - T sin (9j/2) and define the 
scaled modulation function 

f N (0) = F N (T sin 2 (0/2)) = (-iy , (50) 

for 6 e [jn/(N + 1), (j + l)n/(N + 1)]. With 

9 & 

sin 2 " - sin 9 = (20" 2p £ C' 2p sin [0 - r + 1) 8] , (51) 

we can write F PuP2t ... >Pn as a linear combination of terms in the 
form 

U,~ * = \ de "'\ i dQl J ' dGl fl /w W Sin M>) ' 

(52) 

with Lj < pj + 1, Suffices it to show fq l ,q z ,- ,q„ - f° r odd n 
and 2/ = i W/ ^ AT. We notice that fy (9) is a periodic function 
with a period of 2n/(N +1) and therefore can be expanded 
into Fourier series 



/*(#)= J] r sia [k (N +l)6] - 



fc=l,3,5,- 



(53) 



The key feature of the Fourier expansion to be exploited is 
that it contains only odd harmonics of sin[(JV + 1)6], With the 
Fourier expansion, we just need to show that 

dO n --- d9 2 dGi fl cos { r fij + VJ G j) = °' ( 54 ) 
o Jo Jo l j= l 

for n being odd, rj being an odd multiple of (N + 1), and 
2Z%i \lj\ ^ N. With the product-to-sum trigonometric for- 
mula repeatedly used, it can be shown by induction that af- 
ter an even number of variables 9i , 62, . . . , 62k have been inte- 
grated over, the resultant integrand as a function of 02k+i is the 
sum of cosine functions of the form 



COS (R.2k+\62k+l + Q2k+\02k+l) , 



(55) 



with /?2t+i being an odd multiple of (N + 1) and IQifc+il ^ 



2/t+l 



Yj~; = i \qj\- In particular, the last step is 



f 

Jo 



cos(R„9„ + QnO„)d0 n . 



(56) 



Since R„ is an odd (non-zero, of course) multiple of (N + 1), 
and|g«l < Z'Li ?y ^ N, we have R„ + Q„ + and the integral 
above must be zero. Thus Eq. d49l holds. The proof is done. 

It should be noted that in the proof above, the perturbation- 
theoretical expansion requires that the Hamiltonian of the bath 
have a finite norm, which means that the noise spectrum felt 
by the qubit has a hard cutoff. 



4. Universality of UDD against population relaxation 



5. Time-dependent Hamiltonians 



A straightforward corollary of Eq. (|49b is that UDD can also 
be used to suppress population relaxation of the qubit. Con- 
sidering the most general qubit-bath Hamiltonian in Eq. (fTTT i 
and assuming that the UDD sequence consists of N instanta- 
neous ti pulses to rotate the qubit around the z-axis, we aim to 
show that the relaxation of the qubit population in |f ) and |X) 
is eliminated up to 0(T N ). The propagator of the qubit-bath 
evolution from to T is 

U(T,0) =U Q (T - T N )& z O (T N - 7V0 

■■■& z U Q (T 2 -T x )fr z U Q (Ti), (57) 

where Uo(i) - e~ lHt is the free propagator. U(T, 0) can be 
written as Un (for N being odd) or & z Un (for N being even) 
with 

N _ e -i[C'+{-lfbW-T N ) . . . e -«(C'-D)(7' 2 -7' l ) e -/(C'+D)7' 1 ^ ^ 

in which the Hamiltonian has been separated into C' = C + 
fr z ® Z and D = fr x ® X + fr y ® Y. With the definition D, (?) = 
e De , the propagator can be formally expressed as 

U N = e- itT 3-e-'' £ F ^o im = tft**a + fljW (59) 

where 

y(even) 



u 



,v 



--'"X(-/f4 



/,=o 



consists of even powers of D, and 



A=0 



(60) 



(61) 



consists of odd powers of D, with A,', obtained from A„ in 
Eq. (144-b by replacing Zi(t) by Dj(t). By expanding Di{t) into 
Taylor series [similar to Eq. d42b l 



D I (t) = 2 J D P t p , 

p=0 



(62) 



the identity Eq. (ggj) immediately gives A' 2k+l = 0(T N+l ). As 
a result, t/]£ = 0{T N+l ) and the propagator 

C/jv = U%™*> + 0(T N+1 ), (63) 

contains only even powers of D up to 0{T N ). Since D contains 
only the Pauli matrices & x and <x v and an even power of the 
two Pauli matrices d^dy (with n x + «,. being even) is either 
unity or i& z , the propagator 



From the procedures following Eqs. d42b and d62l ). it is 
immediately observed that the proof above applies to time- 
dependent Hamiltonian as long as a Taylor expansion of the 
Hamiltonian similar to those in Eqs. d42l and d62b exists (such 
as a Hamiltonian having analytical time-dependence). Such a 
generalization was presented by Pasini and Uhrig [55]. 



6. UDD with non-instantaneous pulses 

With the help of Eq. d54l i, we realize that Eq. d49t holds 
for more general modulation functions F N (t) as long as the 
scaled modulation function /w(0) = Fn \T sin (6/2)) contains 
only odd harmonics of sin[(A^ + 1)9] as in Eq. d53l . i.e, 

oo 

f N (e) = J]A k sm[(2k+l)(N+l)8], (65) 

with arbitrary coefficients Aj. Motivated by this observation, 
we try to generalize UDD to the case of n pulses with a finite 
duration. 

For the case of UDD against general decoherence, we con- 
sider the control of the qubit by an arbitrary time-dependent 
magnetic field B(t) applied along a certain direction to protect 
the qubit coherence along this axis. Without loss of general- 
ity, we take this direction as the z-axis. The general qubit-bath 
Hamiltonian under DD control is 

H(t) = C + fr x ®£+ fry ®Y + fr z ®Z+ -fr z B(t). (66) 

In the rotating reference frame following the qubit precession 
under the magnetic field, the Hamiltonian becomes 



H R (t) = C + cos[(p(t)]D + + sin[<p(t)]D~ 



(67) 



where the precession angle (p(t) - J B{t')dt' , C = C+fr z ®Z, 
D + = fr x ®X+fr x ®Y, and£>~ = fr x ®Y-fr x ®X. The propagator 
in the rotating reference frame is 

U N = e- iC ' T ,? exp \-i f YF*(t)DfWt\, (68) 



with F + N {t) = cos[<f>(t)], F N (t) = sin[0(f)], and D*(f) 



(even) 



e' c 'D e . Again we decompose Un as the sum of U N 
(which consists of even powers of D*) and U° (which con 
sists of odd powers of D ± ), 



(J - e -iH cB (T)T+0(T N+1 ) 



(64) 



where the effective Hamiltonian H e s(T) - C e g(T) + fr z ® 
Z e ff(T) contains only pure dephasing term fr z ® Z e g{T) and 
commutes with fr z . Thus the A^-pulse UDD eliminates the 
population relaxation up to 0(T N ). 



U 



(even) 
N 



U 



(odd) 



k=0 A.i~A n 



(69a) 



'"I, Z (-O 2 ^ 1 ^^ , (69b) 
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FIG. 2: (Extracted from Ref. 15011 . Copyright of the American Phys- 
ical Society 2008) An example of (a) the scaled modulation func- 
tions f^j(9) for the generalized 3rd order UDD control and (b) the 
corresponding magnetic field B(t). The dashed lines indicate the cor- 
respondence between the sudden jumps of the modulation function 
f^{0) in (a) and the sharp spikes as instantaneous w-pulses in (b). 



where 



= F^(t n )dt n F^(t n -i)dt n - 
Jo Jo 

f 

Jo 



F%(h)dti D A { (t„) Df-< {t n .i) ■■■D*' (h), 



(70) 



initial conditions /^(0) = 1 and f N (Q) - 0, one can choose 
the field such that f N (8) is continuous while /^(0) has sudden 
jumps between +1 and -1 at 8 — jn/(N +1). Thus, a gener- 
alized UDD sequence can be chosen the following way: For 
8 € [0, n/(2N + 2)], fff(8) can be arbitrary but sudden jumps 
from -1 to +1 at 8 - and from +1 to -1 at ;r/(2JV + 2), 
and f N {8) is determined from the normalization condition as 

f N (8) = + yi - [f^(0)f. At other regions, /*(#) are deter- 
mined by the symmetry requirements. The pulse amplitude 
B(t) for the generalized UDD is 



B(t) 



1 



Kit) dt 



F N (t) = Yx8(t- Tj) + BextraC), (72) 



./=' 



which is a superposition of the instantaneous UDD pulses and 
an extra component B emiL (t) being arbitrary but subject to the 
symmetry requirements. The demand of <5-pulses in the gener- 
alized UDD is consistent with the previous finding in Ref. Il56ll 
that the effect of an instantaneous 7r-pulse on the evolution of 
a qubit coupled to a bath cannot be exactly reproduced by a 
pulse with a finite magnitude. An example of the scaled mod- 
ulation functions and the corresponding magnetic field for the 
generalized 3rd order UDD control are shown in Fig. [2] No- 
tice that due to the variable transformation from 8 to t, the 
magnetic field B{t) does not have the symmetries as the scaled 
modulation functions ft(0). For example, B(t) is not periodic 
and the pulse at different time has different width. 

Obviously, the same argument holds for DD against pure 
dephasing just by changing the rotation axis. 



has a structure similar to A„ in Eq. d44l . After expanding 
[DMtk)) into Taylor series, the identity Eq. d49l immediately 
gives AL'.'j 2 = 0(T N+1 ). Thus the qubit decoherence along 
the z-axis is suppressed to 0\T N+1 J] as long as the scaled 

modulation function fi(8) = F* \T sin 2 (6'/2)J contains only 
odd harmonics of sin[(JV + 1)8] as depicted in Eq. d65l l. This 
condition is satisfied if and only if the scaled modulation func- 
tions fj!j(ff) have the following symmetries: 

1. periodic with period of 2n/(N +1); 

2. anti-symmetric with respect to 9 - jn/(N +1); 

3. symmetric with respect to 8 - (j + l/2)n/(N + 1). 

The anti-symmetry condition requires fi(0) be either zero 
or discontinuous at 6 = jn/(N + 1). But f^(8) and f^(8) 
cannot be simultaneously zero since they have to satisfy the 
normalization condition 



lf^(0)f + UnW] 2 = 1, 



(71) 



according to the definition of F^(t). So there must be sud- 
den jumps at least in one of two modulation functions at 
8 = jn/(N +1), which means the controlling magnetic field 
B(t) has to contain a 5-pulse for ^--rotation at t — Tj. With the 



7. UDD with pulses of finite amplitude 

In realistic experiments, the pulses have finite durations and 
amplitudes, which introduces additional errors. There is a no- 
go theorem which states that instantaneous 7r-pulses cannot be 
approximated by pulses of finite amplitude and of short dura- 
tion t p with accuracy higher than the order 0(t p ) without per- 
turbing the bath evolution 15a , 15711 . However, as we have dis- 
cussed above, the symmetric requirements of fi?(8) automat- 
ically guarantee the performance of UDD. Uhrig and Pasini 
showed that by appropriately designing the pulses, the qubit- 
bath Hamiltonian describing pure dephasing can be trans- 
formed into the form 15811 



H = C + F N {t)cr z ® Z + 0(O, 



(73) 



with the modulation function taking values from {-1,0, 1). 
The scaled modulation function /n(8) = F^ [T sin 2 (#/2)J 
is designed to have the symmetries required in the previous 
proof and therefore can be expanded by odd harmonics of 
sin[(Af + 1)8]. Thus, the decoherence is suppressed up to the 
order 0{T N+1 ) + O(t^). This sequence can also suppress lon- 
gituddinal relaxation Il50l 15811 . An arbitrary order M of pulse 
shaping can be achieved by a recursive scheme based on con- 
catenation [591. 
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D. Comparison of decoupling efficiencies of UDD and CDD 

We consider a DD sequence of N pulses, with a total 
evolution time T and a minimum pulse interval t. In CDD, 
the decoupling order is « ~ log 2 JV and t = 772". In UDD, 
the decoupling order is N and r ~ T/N 2 . To be specific, 
our discussion is based on the pure dephasing model. The 
situation for the general decoherence model is similar. We 
compare the efficiencies of UDD and CDD in the two 
following scenarios: 

Case I: The total evolution time T is fixed. The decou- 
pling precision (defined as the effective coupling under the 
DD control relative to the original one) in UDD was derived 
aslHl 



euDD ~ \\H\\ N T N /N\. 



(74) 



In CDD, it scales with the time and the decoupling order 
as HHH 



«CDD~(ll#l|77Viv) n = (||i/||7T/2 



n 2 /2 



(75) 



Thus with T fixed, increasing the decoupling order and 
hence the number of pulses always increases the decoupling 
precision. An arbitrarily high decoupling precision can be 
achieved simply by choosing a sufficiently high order of DD 
(and correspondingly a sufficiently small pulse interval r). In 
the high-fidelity regime (T is small), the decoupling precision 
of UDD scales with the number of pulses much faster than 
that of CDD. However, if we compare the efficiency of UDD 
and CDD of the same decoupling order n, i.e., the rath order 
UDD (containing n pulses) and the nth order CDD (contain- 
ing 2" pulses), CDD has a much higher decoupling precision 
than UDD does (T"/2" 2/2 <z T"/n\ for large «), since the 
minimum pulse interval t = T/2" in CDD is much smaller 
than that in UDD (r ~ T/n 2 ). For the same reason (namely, 
reduction of t), to achieve a given order of precision, CDD 
indeed requires by far less than the seemingly exponential 
cost. 

Case II: The minimum pulse interval r is fixed, which 
is a frequently encountered restriction in realistic experi- 
ments. In this situation, increa sing the order of DD leads 
to two competing effects 11521 16011 . First, the qubit-bath 
coupling is eliminated to a higher order, which tends to 
increase the decoupling precision. Second, the total evolution 
time T increases and the bath has more time to inflict qubit 
decoherence. Competition between these two effects leads to 
the existence of an optimal decoupling order, beyond which 
further increasing the order of DD does not improve the 
decoupling precision any longer. For a given minimum pulse 
interval r, the optimal order of UDD is IJ52I1 



«opt,UDD~l/(l|tf||)T), 

and that of CDD is (MSI H 



«o P t,CDD ~ -log 2 (||//||r). 



(76) 



(77) 



UDD has a much higher optimal level than CDD for a small 
minimum pulse interval, and therefore the highest decoupling 
precision that can be achieved by UDD is much higher than 
that by CDD. 



E. Experimental progresses 

UDD was first realized in experiments by Biercuk et al. in 
an array of ~ 1000 Be + ions in a Penning ion trap J40l |6U 16211 
with noises mimicked by artificially introduced random mod- 
ulation of the control fields. The qubit states were realized 
using a ground-state electron-spin-flip transition. Coherent 
qubit operations were achieved through a quasi-optical mi- 
crowave system. UDD was compared with CPMG in the "low 
fidelity" regime for various classical noise spectra. The data 
show that UDD dramatically outperforms CPMG for Ohmic 
noise [power spectrum S (u>) oc a;] with a sharp cutoff, while 
for the ambient magnetic field fluctuations whose power spec- 
trum S (a>) oc 1 /(jj 1 has a soft cutoff, UDD performs similarly 
to CPMG over the entire range of accessible noise intensities, 
consistent with the theoretical predictions [38, 53]. 

The first experimental realization of UDD against realis- 
tic noises was achieved by Du et al in a solid-state system, 
namely, irradiated malonic acid single crystals. The spins of 
the radicals in the crystals created by irradiation form an en- 
semble of qubits. The nuclear spins, in samples with rela- 
tively low concentrations of radicals, constitute the quantum 
bath which can be considered as finite for the time-scales in- 
volved in the experiment and therefore has a finite noise spec- 
trum. The pulsed electron paramagnetic resonance was used 
to demonstrate the performance of UDD for preserving elec- 
tron spin coherence at temperatures from 50K to room tem- 
perature 14111 - Using a seven-pulse UDD sequence, the elec- 
tron spin coherence time was prolonged from 0.04 ms to about 
30 ms. The experimental data from different samples un- 
der various conditions demonstrate that UDD performs bet- 
ter than CPMG in fighting against noises from nuclear spins. 
The good agreement between the experiment and calculations 
based on microscopic theories llm \Y% enables the authors to 
identify the relevant electron spin decoherence mechanisms 
as the electron-nuclear contact hyperfine interaction and the 
electron-electron dipolar interaction. 



III. NEW DEVELOPMENTS 

A. CUDD: Concatenation of UDD 

CDD can eliminate all the qubit-bath couplings (including 
pure dephasing and population relaxation) up to an arbitrary 
order N at the cost of exponentially increasing number (of the 
order 4") of controlling pulses. In contrast, UDD sequence 
uses the least number (i.e., N) of controlling pulses to elimi- 
nate either pure dephasing or population relaxation (but not 
both) to the desired order N. Based on a combination of 
CDD and UDD, a new DD sequence (named CUDD) was pro- 
posed [63] to suppress both the pure dephasing and the pop- 



10 



ulation relaxation to order N with a much less (of the order 
M^) number of pulses. The essential idea of CUDD is to use 
the Mh order UDD sequence (instead of the free evolution) as 
the building block of CDD sequence. 

The propagator £//v-udd(70 for the qubit-bath evolution 
driven by the general Hamiltonian Eq. (TTTb under Mh order 
UDD sequence of n rotation around the z axis is 

JVudd = e-^'-K-D^KT-TV) . _ e . K(: --bXT 1 -T i ) e -KC' + 6)T i 



-iH cB (T)T+0(T N+l ) 



(78) 



[see Eq. J58H . where {Tj} are given by Eq. ( fT9b and H e ff{T) - 
C e fi(T) + fr z ® Z e ff(r) is a pure dephasing Hamiltonian. The 
pure dephasing can be eliminated by embedding Un-\jdd into 
the structure [& x {- ■ ■ )& x ] (•••)■ The propagator for the with 
order concatenation of Un-vdd is 



V 



(m) 
N-UDD 



<Tx U /V-UDD "-*' U N-UDD ' 



(79) 



with £/^! UDD = t/w-uDD- In the CUDD scheme, £/^° UDD 
eliminates both the pure dephasing and population relaxation 
up to the Mh order with 0{N2 N ) pulses. 



B. Near optimal DD by nesting UDD 

Recently West et al proposed a near optimal DD (64 
obtained by nesting UDD sequences, dubbed quadratic DD 
(QDD), to protect qubits against general noises. The inner 
Mh order UDD eliminates population relaxation and the outer 
Mh order UDD eliminates the pure dephasing, so that both 
pure dephasing and population relaxation are eliminated up 
to the Mh order of the evolution time. Using U n _ udd {t) to 
denote the qubit-bath propagator driven by the general Hamil- 
tonian Eq. ( TTTb under the Mh order UDD sequence of n ro- 
tation around the z axis, the propagator of the (N, M)th order 
near optimal DD, 



U 



(M) 



--U 



(Z) 



>(Z) 



M-\ 



n-udd^T Tm)o- x U n _ vdd (Tm 
■vM'(T 2 - W.tf® (TO 



(80) 



is obtained from Eq. ( IBTl i by replacing the free propagator 
U (t) by U ( ?! VDD (t), where {7^} are given by Eq. <[T9]) with N 
replaced with M. Thus U N eliminates both the pure dephas- 
ing and the population relaxation up to the Mh order using 
0{N 2 ) pulses. Numerical simulation shows that for a fixed 
number of pulses, this DD sequence outperforms CDD and 
CUDD by exponential saving of the number of the pulses and 
it is nearly optimal for small N, differing from the optimal 
solutions by no more than two pulses. 

A proof of the QDD was attempted in Ref. 15511 with the ar- 
gument that after the inner level of UDD control, the resulting 
effective Hamiltonian is time-dependent and the outer level of 
UDD control applies to time-dependent Hamiltonians. The 
effective Hamiltonian under the inner level of UDD control as 
defined in Ref. [55], however, is only piecewise analytical. It 
can be shown by some counter examples 116511 that for a gen- 
eral piecewise analytical Hamiltonian taken as resulting from 



certain inner level of control, it is not guaranteed that the outer 
level of decoupling can be realized to the desired order. Thus 
it remains an open question to us why the nested UDD control 
works. 



C. Protecting multi-qubit states by UDD 



Mukhtar et al recently showed H66H that by applying a se- 
quence of unitary operations 



P* = 2\m#\ ~ A 



(81) 



on the multi-level quantum system according to the timing of 
UDD, the initial quantum state \ip) is protected to the order of 
O \T N+[ ). This operation was also given in Ref. 15111 . 
Obviously, we have P'. - P^. We define the operators 



C = (H + P^HP^a, 
Z = (H - P^HP f )/2. 



(82a) 
(82b) 



Then the system-bath Hamiltonian is separated into two parts 
H = C+Z, (83) 

where C commutes with the operator P,/, while Z anti- 
commutes with P.,., i. e, 



•i" 



Pij/CPij, — C, 
Pif/ZP,/, — —Z. 



(84a) 
(84b) 



By applying a sequence of Af operations P^ according to the 
timing of UDD, the system-bath propagator reads 



U N =p N e - iiC+z)(T - TN} P ll ,e~ KC+z)(TN ~ TN - ,) 

...p^C^)(T2-T0p^-i(C + Z)T^ 



(85) 



Note that a final P^ pulse is required for odd M Similar to the 
procedure in the proof of the universality of UDD, we rewrite 
the propagator as 



U N = e -iCT£ e -i£F N (t)Z,( f )dt > 



(86) 



where Z/(f) = e' c 'Ze ,Cl anti-commutes with P^,. We separate 
On into two parts 



where 



fj _ /V(even) /V(odd) 

Un = U N + U N , 

CO 
CO 

f, (odd) _ -iCT V ( ;\ 2k A 

U N -e 2-T ' 2t+1 ' 



(87) 

(88a) 
(88b) 



/(=o 
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with 



A„ = I F N {t„)dt„ I F N (t„-i)dt„-i ■ ■ ■ 
Jo Jo 

F N (h)dhti(t H )Z I (t„-l)---Z I (h). 

Jo 



(89) 



Obviously, U^ commutes with Py, since it contains even 
powers of Z. Following the same arguments in the proof of 
UDD for qubit dephasing, we conclude that U° - 0(T N+l ). 
Thus, P^Un - UnPip + 0(T N+l ), which immediately indicate 
that the expectation value of P^ and hence the quantum state 
\ijj) are preserved up to 0(T N+l ). 



IV. CONCLUSIONS AND PERSPECTIVES 

In summary, we have given a review of recent progresses 
in protecting qubit coherence by the dynamical decoupling 
schemes. The DD techniques are originated from the mag- 
netic resonance spectroscopy. The developments for quan- 
tum information technologies can in turn advance the high- 
precision magnetic resonance spectroscopy. For example, 
UDD has recently been applied in magnetic resonance imag- 
ing of tumors in animals 16711 . Extension of the spin coherence 
by DD may have important applications in nano-scale or even 
atomic scale magnetometry 16811 . 

Remarkably, experiments have demonstrated the DD 
method as a particularly promising scheme for protecting 
quantum coherence in quantum computing. As compared to 
the quantum error correction schemes, the DD requires no 
auxiliary qubits and can be integrated naturally with the quan- 
tum gates without extra hardware overhead. However, the DD 
approach has a shortcoming in that it works only for slow 
baths or for non-Markovian noises, in the sense that the char- 
acteristic separation time of the DD sequence is required to 
be shorter than or at least comparable to the inverse of the 
characteristic width of the noise spectrum. The quantum er- 
ror correction scheme has no such requirements. In dealing 



with errors in quantum computing due to spontaneous emis- 
sion, combination of DD and quantum error correction was 
proposed 16911 . It is conceivable that in future quantum com- 
puting, the non-Markovian noises be decoupled by DD and 
the remaining Markovian noises be coped with by quantum 
error correction. In general, for a multi-qubit system coupled 
to both Markovian and non-Markovian noises, a combination 
of the two paradigmatic error-countering methods provides a 
complete picture for scalable quantum computing 114411 . 

In the present research of DD, mostly the pulses are as- 
sumed instantaneous with only a few exceptions. Two im- 
portant issues are under intensive research, and some remark- 
able results have emerged recently JH H2, [3, HH HI] . One 
is how to extend the DD to implement high-fidelity quantum 
gates or hybrid DD with quantum gates. Can some ideas be 
borrowed from DD for realizing dynamical control resilient 
to noises? Such an issue was previously addressed in simu- 
lation of quantum processors with DD approaches 17011 . Re- 
cently, encouraging progresses have been made toward hy- 
bridization of quantum gates and DD J42l - l44ll . Another issue 
is how to design a quantum gate (such as a qubit flip, which 
is required in DD) optimally in the presence of environmental 
noises. Various optimization schemes have been invented for 
suppressing/minimizing the noise effect to a certain order Il7ll - 
17411 . Ref. 15911 has established a systematic method to achieve 
an arbitrary order of precision based on iterative construction 
of finite-amplitude pulses. It is of interest to ask whether and 
how the pulse shaping for quantum gates with an arbitrary 
order of precision can be systematically constructed without 
iteration, with the development from CDD to UDD being an 
inspiring example. 
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